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Abstract. We investigate through Monte Carlo simulations the non-equilibrium behaviour of the three- 
dimensional XY-model quenched from a high temperature state to its ferromagnetic and critical phases. 
The two-times autocorrelation and response functions are determined in the asymptotic (scaling) regime, 
from which the nonequilibrium exponents A and critical A c are extracted. The form of the scaling function 
is in agreement with the prediction of local scale-invariance. The so-called limit fluctuation-dissipation 
ratio Xao is shown to vanish in the ordered phase and to reach a constant value around 0.43 for the critical 
quench. 

PACS. 75.40.Gb Dynamic properties 

05.70.Ln Non-equilibrium and irreversible thermodynamics 
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1 Introduction 

After the investigations on aging in glass models. [1|2I3| 
efforts have been concentrated on the study of non- 
disordered model systems since it was soon realized that 
some of the basic aspects of aging are present in sim- 
ple systems too with some caracteristic features. |4I5| In 
most studies, special attention is put on pure ferromag- 
netic models undergoing a second order phase transition. 
Typically, one cools down a sample initially prepared in 
its high temperature phase to its ordered phase. As time is 
running, domains start to grow with a linear size I ~ t x l z . 
In the thermodynamic limit, equilibrium is never reached 
which implies an infinite relaxation time.|S] In order to 
probe aging, which implies a full dependence on the past 
evolution, it is useful to compute two-times correlations: 



C(t, t u 
and linear responses 

R(t) t w 



(a(t)a(t w )) 
S(a(t)) 



Sh(t v 



h=0 



where a is the order parameter, t the observation time and 
t w the waiting time (< t). 

Usually, at a coarse-grained level, the dynamics of such 
systems is described by a Langevin equation: [Hj 



dt 



SH\ 



6<f>i 



■Vi(t) 



(3) 



where H is the free-energy functional, f]i(t) is a thermal 
gaussian noise at site i with: 



M«)> = , (vi(t)v j (t')) = 2TS(t-t')6 i , j 



(4) 



As mentioned above, the breakdown of time translation in- 
variance implies that the two-time functions do not merely 
depend on t — t w but explicitly on t and t w . On the bases 
of general scaling arguments, one postulates the follow- 
ing scaling laws for the autocorrelation and autoresponse 
functions: 

C(t,t w )^t- b f c (t/t w ) (5) 



and 



R(t,t w )fat w a fR.{t/t w ). 



(6) 



In the asymptotic limit, where t 3> t w ^> 1, fc an d Jr 
decay algebraicly: 



fc(x) 



-\ c /z 



Sr{x) 



-Ah/* 



(7) 



where z is the dynamic exponent and Xc,r are respectively 
the correlation [7] and response exponents. |33| 

A consequence of the breakdown of time-translation 
invariance is the violation of the fluctuation-dissipation 
theorem (FDT) which relates for equilibrium systems the 
response to the correlation function via: 



R(t — t w ) 



1 dC(t - t w ) 



T 



(8) 



When the system is not at equilibrium R{t, t w ) and 
C(t,t w ) both depend on t and t w , and the FDT JSJl does 
not hold anymore. One introduces a new function X(t, t w ) 
defined by: [H] 



R(t, t v 



x(t,t w ) dc(t,u 
t dt w 



(9) 



2 



S. Abriet and D. Karevski: Off equilibrium dynamics in 3-d XY system 



which measures the deviation from equilibrium. In some 
mean-field theory of glassy systems, where it was first for- 
mulated, the fluctuation-dissipation ratio (FDR) X(t,t w ) 
depends on time only through a functional dependence 
on C(t, fu,),P] that is X = X(C(t,t w )). In this case, 
integrating with respect to t w , one gets T\ — 
Sc(t t ) X{C')dC' , so that the slope of the parametric plot 
susceptibility \ versus the correlation function C will give 
the FDR. In the asymptotic regime, for such pure ferro- 
magnetic models, the FDR reaches a constant value X^ 
after an initial quasi-equilibrium regime where X = 1 
(FDT valid). It was argued that this asymptotic con- 
stant is universal and characterises the non-equilibrium 
process. |llj When a ferromagnetic system is quenched to- 
ward its ferromagnetic phase, the long-time dynamics is 
governed by the diffusion of domain walls and the coars- 
ening of domains. This coarsening leads to a vanishing 
FDR X x = 0. [5I8I9I10| . The scaling functions defined 
previously are fullfilling this requirement. In the case of 
a critical quench, that is exactly at the critical tempera- 
ture, several studies on various systems have given some 
insights. In that case, the reason for aging is related to 
the development of the spatial correlations over a length 
scale £ that grows as t 1 '*" where z c is the critical dynamic 
exponent. So the system is still disordered over a length 
scale larger than while it looks critical in regions 

smaller than £(t). In the thermodynamic limit, the equili- 
bration is unaccessible since the system is of infinite size so 
£ indefinitely increases. For a fully disordered initial state, 
one has in the kinetic spherical model X^ = 1 — 2/d for 
2 < d < 4 and = 1/ 2 ford > 4,[TT] in the Id Glauber- 
Ising model 1^ = 1/2. [nil's*] For algebraic initial corre- 
lations, it was shown on the kinetic spherical model that 
there exists a rich kinetic phase diagram, depending on 
the space-dimension d and the algebraic decay exponent 
ii>-|33) where the limit value J M is either vanishing or 
depending on d and u). In the Id Glauber-Ising model, 
it was shown that X^ is independent on o;.[13| A sys- 
tematic field-theoretical approach was developped to cal- 
culate the limit FDRjT^ and applied to several models, 
ferromagnetic systems at criticality, dilute Ising model, 
model C dynamics. [151161 IT] Simulations in 2d systems 
gave Xoo = 0.34 for the Ising model. |3~1I18| X^ = 0.41 
for the three states Potts model, Xoo = 0.47 for the 4- 
states Potts model. [22] In 3d, X^ = 0.40 was obtained for 
the Ising model. [HQ Recently, performing a Monte Carlo 
simulations with Glauber dynamics, we found from the 
susceptibility versus correlation function slope a continu- 
ously varying FDR in the 2-d XY model quenched below 
the Kosterlitz-Thouless point. [H"] But in this case, one has 
to use with great care this slope, since the actual defini- 
tion of the FDR X^ leads to a logarithmicaly vanishing 
result. In this paper, we present the results obtained on 
the 3-d XY model quenched onto criticality and below. 



2 The 3d XY-model 

The XY-model is one of the most popular theoretical 
models. As it is well-known, in the two-dimensional case 



there is no finite magnetization at non-zero tempera- 
ture, but nevertheless there is indeed a topological phase 
transition occuring at finite temperature due to vortex 
pairing. 20 21 In three dimensions, the scenario is more 
conventional since the additional dimension permits the 
existence of an ordered phase at finite temperature. Thus, 
the 3d-system presents a second-order phase transition 
which is probably related to the density of vortex-strings 
22 23] . The hamiltonian reads: 



(i,3) 



s, 



(10) 



where i, j are nearest neighbour sites on a cubic lat- 
tice and the Si are two-dimensional classical spins whose 
length is set to unity. From high-temperature series anal- 
yses, the predicted critical inverse temperature is |24| : 
[3 C = 1/T C = 0.4539 ± 0.0013. In our simulations, we use 
the value (3 C = 1/T C = 0.4542 obtained from a Monte 
Carlo study [25] • The critical exponents have been es- 
timated by renormalisation group techniques and e- 
expansions. [26127] In our study we take the values of the 
critical exponents obtained in ref.|28j from a numerical 
approach: 

/3~ 0.349 ^- 0.672. (11) 

Simulations on the three-dimensional XY model with pe- 
riodic boundary conditions have revealed that the critical 
dynamics exponent z c is close to 2 [29130) . For a quench be- 
low the critical temperature, that is in the ordered phase, 
one has z = 2.[f"j 

3 Numerics 

We simulate the dynamics through a single spin-flip algo- 
rithm of Glauber like form where the transition rates p 
related to the single flip {6i} — > {#■}, where the 9s are the 
angular variables of the XY model, reads: 



exp(-/3W[6» l ]) + exp(-/3H[0<] 



(12) 



with (3 the inverse temperature. One may notice here 
that the Metropolis algorithm, with transition rate p = 
min[l, exp(— (3AE)], leads to the same time evolution 
of the autocorrelation and susceptibilities, as we have 
checked it. 

The two-times autocorrelation function is calculated 



via 



C(t, t v 



L' :i 



(13) 



where the brackets ( ■ ) stand for an average over initial 
configurations and realisations of the thermal noise. 

The two-times linear magnetic susceptibility x(t,t w ) 
can be computed by the application of a random magnetic 
field. The amplitude of the field has to be small in order to 
avoid nonlinear effects. The ZFC-susceptibility is usually 
computed utilising a random distributed field via:[8j 



X(t, t-u 



ry-3 y^^K cos 9i(t) + h y sin 9i(t)) 



h 2 L 3 



(14) 
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where the line stands for an average over the field realisa- 
tions. Indeed, here we concentrate on a different approach 
which is based on the lines of ref.|3J. By definition, the 
linear autoresponse function R(t,t w ) is 



R(t, t v 



1 

L3 



E 



5{cos6i(t)) S(sm9 t {t)) 
Sh x (t w ) 6h y (t w ) 



h=0 



Following ref.|31|. we can reexpress the linear response 
function in the form:[T§| 

R(t,t w ) = /3(cos 6i{t) [cos 6i(t w + 1) - cos6?(t w + 1)]) 

+0{am6 i (t)[faRB i (t vl + 1) - sm6?{t w + 1)]> , (16) 

where cosOf and sin#™ are the components of the local 
Weiss magnetisation: 



S 



1 d 



(3 dh 



■ In Zi 



x.y 



(17) 



h=Q 



where Z l = exp (-/3H(6 t , h)) + exp(-f3H(8' i: h)) is the 
local partition function in the field. Finally the fluctuation 
dissipation ratio is calculated via 



X(t, t w ) 



TR(t, t u 



c(t,t w + 1) - c{t,t w ) ■ 



(18) 



4 Results 

4.1 Quenches below T c 

We present the results obtained on systems of linear size 
L = 50, with periodic boundary conditions, quenched in 
the low-temperature phase from an infinite-temperature 
initial state. The averages are performed with the use 
of more than 5000 samples. The expected behaviour of 
the two-times autocorrelation function for coarsening pro- 
cesses is at sufficiently long times [H]: 



C{t,t w ) = M 2 e JT)f c 



(19) 



where M eq (T) is the equilibrium magnetisation at the 
temperature T. In order to check this behaviour, we have 
simulated quenches at two different temperatures, T = 0.9 
and T — 1.5, using the Glauber like dynamics. The results 
are shown respectively on figure 1 and 2. 

For small time separations, we clearly see that the 
correlations rapidly decay to the expected plateau values 
M^ q (T) , here estimated numerically by the use of an inde- 
pendent Wolff algorithm simulation. |34I35| For large time 
separations, the system enters in a scaling region where 
the correlation functions are described by a scaling func- 
tion fc of the variable x — t/t w . As expected, the aging 
scaling function finally decays algebraicly in the asymp- 
totic limit with an exponent Ac roughly about 1.7 ± 0.1. 
However, we clearly see that as the temperature is in- 
creased the finite-size effects become stronger. 
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Fig. 1. Asymptotic behaviour of the autocorrelation function 
at T = 0.9. The dashed line corresponds to the value M^ q (T) 
and the solid line represents the algebraic decay of f(t/t w ). 
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Fig. 2. Asymptotic behaviour of the autocorrelation function 
at T = 1.5. The dashed line corresponds to the value M^ q (T) 
and the solid line describes the algebraic decay of f(t/t w ). 



If one accepts the value Ac = 1.7, since it is very close 
to the free field value A c = d/2 = 1.5, one can expect a 
scaling form of the two-times autocorrelation function not 
too far from the free- field form:|3Hj 



C°(t,t w ) = M? q (T) 



2 \ --^c/ 2 



t/t v 



(20) 

Accordingly, we have replotted in the inset of figure 2 the 
two-times autocorrelation as a function of the scaling vari- 
able y = . Within this new scaling variable, the 
power-law form, with Ac = 1.7, is remarkable. 

In the low-temperature phase, the behaviour of the 
two-times response function had been conjectured to be 
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(see Ref. 0] for a review): 

i?(M™) ^ C~ a /* (21) 

where the scaling function fn{x) rs x~ Xr / z at x ^> 1. 
Furthermore, if one assumes that the response function 
transforms covariantly under local scale-transformations, 
one expects for the scaling function fa |37I38| 

f R (x) =r x 1+a - XR / z (x-l)- 1 - a (22) 

where ro is a normalization constant. For a disordered ini- 
tial state we expect \r = Ac = A, with the nonequi- 
librium exponent A bounded bv |39l7l40| : d/2 < A < d, 
where d is the euclidian dimension of the system. 

We have computed the autoresponse function at tem- 
peratures T — 0.9 and T = 1.5. On figure 3 and fig- 
ure 4(left), assuming that a — 1/z — 1/2, 0T] the rescaled 
response function is plotted as a function of t/t w . The 




10 - j I , , — , — , , — , — 

1 10 100 

t/t w 

Fig. 3. Scaling behaviour of the response function at T = 0.9 
for different waiting times. The solid line gives the local-scale 
invariance prediction with ro = 0.18. The dashed line is a guide 
to the eyes for the power law behaviour. 

collapse of the different waiting time curves is quite good. 
The value a = 1/2 had been already found numerically 
and analytically in the 3D Glauber-Ising model 38 9 and 
the 3D kinetic spherical model |33lllj . We have performed 
an analysis at fixed ratio t/t w that corroborates the power 
law prefactor in (I21|l with a close to 1/2 at large t/t w as 
seen on figure 4. In fact, numericaly we obtain a value 
which is slightly larger than 1/2, around 0.6 but the ten- 
dency when increasing t/t w is toward the decrease of the 
effective a value. 

Due to important thermal fluctuations the algebraic 
regime x~ Xr / z is hardly accessible. However, our datas 
seem to validate the relation Ac = Xr as it can be seen on 
figure 3 and figure 4 for T = 0.9 and T = 1.5 respectively. 
Moreover, as it can be seen from figure 3 and 4, the scaling 




Fig. 4. Left: Scaling behaviour of the response function at 
T = 1.5 for different waiting times. The solid line gives the 
local-scale invariance prediction with ro = 0.18. The dashed 
line is a guide to the eyes for the power law behaviour. Right: 
Power-law decay of the response function at fixed x — t/tw 
ratios and T — 1.5. The results are corroborating the value 
a = 1/2. 

form 1)22(1 predicted from local scale invariance is perfectly 
reproduced by our datas with ro = 0.18. 




20 40 60 



t/t„ 

Fig. 5. Rescaled Fluctuation-Dissipation ratio Xt^ for a 
quench in the ferromagnetic phase at T = 0.9. The inset is 
a parametric plot of X versus correlation function. 

As discussed in the previous section, the asymptotic 
value of the fluctuation-dissipation ratio Xac is expected 
to vanish in the low temperature phase. Assuming that 
Ac = Afl, the fluctuation-dissipation ratio evaluated from 
the asymptotic behaviour of ((19JI and (|21(l vanishes as: 

X(t,t w )~t- a f x (t/t w ) (23) 

where fx{x) is a scaling function related to fc{x) and 
fn(x) with lim-r^oo fx(x) — const., a time-independent 
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constant. On figure 5 we have plotted the rescaled FDR 
obtained numerically at T — 0.9, assuming a = 1/2, for 
different waiting times as a function of t/t w . We see a 
very good collapse of the different waiting time datas. The 
solid line is the expected scaling form obtained from the 
definition of the FDR together with the scaling functions 
fc and fu assuming the local scale prediction with Tq = 
0.18. Asymptotically, the scaling function fx reaches a 
constant value: 



le+01 



lim fx(x) = 



Trn 



(24) 



which is close to 0.08 at T = 0.9. In the inset graph of 
figure 5, we have represented the FDR X(t, t w ) as a func- 
tion of the autocorrelation function C(t,t w ) for different 
waiting times, where we see the vanishing of the X ratio 
as the correlation decay. 



4.2 Critical quench 



u 



o.oi 



0.001 




Fig. 6. Rescaled autocorrelation functions for a critical quench 
of the 3D XY-model. The slope of the solid line stands for 
-\ c c/z c ~ -1.34. 



We focus now on the critical quench starting from a 
high-temperature initial state. In the vicinity of the crit- 
ical point, the magnetisation at equilibrium behaves as 

|T — Tel 13 ~ ^e/^ U and with £ eq ~ t L J ] z " , one expects for 
the two-times autocorrelation function the form: 



C(t,t w ) = A c c t- 2 ^fz(± 



(25) 



The expected scaling for the autoresponse function is 
given by 



R(t, t u 



Ar t. 



-1-2/3/vZc fc 



t 



fn(- 



(26) 



The scaling functions /§ R are expected to decay alge- 
braicly at large time separation with the same exponent 
X C R = A£., however different from the low-temperature one. 

The numerical calculations have been performed on 
periodic cubic lattices of linear size up to L — 50 and the 
physical quantities have been averaged over 5000 noise 
realisations. On figure 6 we have plotted the rescaled two- 
times autocorrelation functions for different waiting times. 
The collapse of the rescaled functions is very good, in 
agreement with the form l|25l) . Moreover, at large t/t w 
we clearly see a power-law behaviour from which by an 
algebraic fit we extracted the exponent X c c j ' z c ~ 1.34 that 
fulfills the scaling bounds d/2 < X c c < d (remember that 
z c is close to 2). Analysing more carefully the datas, we 
obtain a serie of estimates for the exponent Xq/z c that 
are represented on figure 7. Our final estimate is 



X c c /z c = 1.34 ±0.05 



(27) 



We have also explored the finite-size effects in more de- 
tails. For that purpose we have studied the time evolution 
of cubic systems of linear sizes L = 10, 20, 30 and 40. As- 
suming that 



C{t,t w ,L) = b- 2f *l v C{b-*°t,b- Zc t w ,b- x L) 



(28) 



u .4. .8. . s* =i= = 
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t: 



2 



Fig. 7. Effective exponent \ c c / z c for different waiting times as 
a function of t obtained from systems of size L = 50. 



with b — L we obtain the scaling form 
C(t,t w ,)=L 2f3/v T f— , 



(29) 



On figure 8, we have plotted J- for different sizes at fixed 
ratio x = t w /L Za — 0.1 as a function of t/L Zc . Taking into 
account (|25|l and l|29|l , we expect that the scaling function 
T has an algebraic decay with the exponent —\ c c jz c for 
t not too large. The collapse of the curves is quite good. 
For earliest times the system is still in quasi-equilibrium 
and for biggest times the growth of correlated domains 
is limited by the size of the system. The numerical datas 
seem to validate the scaling assumption (|29|) . 

We focus now on the autoresponse functions computed 
with formula i|16|) . The quantities are quite noisy since 
thermal fluctuations are of the same order as the ampli- 
tude of the response at long times. From ©, one should 
have a = b and = = A c since the initial state is 
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l=3o : 
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y~x _ 











t/L' 



Fig. 8. Rescaled autocorrelation functions for different system 
sizes. All the curves collapse in the full-aging regime where the 
scaling function F decays as (t/L Zc ) x c /Zc 



fully disordered. On figure 9, the rescaled autoresponse 
is represented as a function of t/t w for several waiting 
times t w . We see that the decay of the autoresponse func- 
tion is compatible with the algebraic assumption with 
X c R /z c = X c /z c ~ 1.34. 



10 



id- 



io 
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have plotted the FDR X(t,t w ) as a function of t w /t cal- 
culated on a system of linear size L = 30, and for different 
waiting times. In the long time limit t w /t 0, the FDR 
approaches a constant value, around 0.4, but due to the 
numerical noise, basicaly related to the numerical deriva- 
tive of C(i, t w ), it is difficult to give a more precise value. 
However, if one supposes an asymptotic linear form 



X(t,t w ) ~ X 
it is possible to extract the value: 



T ' 



X- 



0.43 ±0.04 



(30) 



(31) 



In ref. the limit FDR was calculated in an e ex- 

pansion. The two-loop expansion leads to a value = 
0.416(8) that perfectly fits, within the error bars, our es- 
timate. 




Fig. 10. Fluctuation Dissipation ratio at the critical point 
as a function of t w /t for L = 20 (open symbols) and L = 30 
(filled symbols) . The solid line gives the analytical result Xoo — 
0.416(8) of ref. |15| . The inset shows the dependence of X with 
the same waiting times as a function of the correlation C. 
e = A-d 



Fig. 9. Rescaled autoresponses as a function of t/t w for dif- 
ferent waiting times and for a critical quench. The solid line 
gives the scaling function as predicted by local scale invari- 
ance. The dashed line gives the algebraic decay (t/t w )~ x R^ Zc 
with X c R /z c = \ c c /z c = 1.34. 

As stated previously, the out-of equilibrium behaviour 
is somehow characterised by the FDR X(t,t w ), calcu- 
lated numerically via l|18|) ■ Assuming the scaling forms of 
the two-times autocorrelation and response functions, the 
FDR is supposed to be a function of the ratio t/t w only. 
This behaviour is well reproduced on figure 10, where we 



5 Conclusion and outlooks 

We have studied the relaxation behaviour of the 3D XY- 
model quenched from its high temperature state to its 
ferromagnetic phase and at criticality. Through Monte 
Carlo simulations, we have computed the autocorrelation 
function and its conjugate response function on cubic lat- 
tices of maximal linear size L = 50. From our datas, we 
have confirmed the general scaling scenario at and be- 
low criticality. The autocorrelation non-equilibrium expo- 
nent Ac below T c was found to be around 1.7. The scal- 
ing analysis of the linear autoresponse function confirmed 
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the equality \r — Xc — A, expected for quenches from a 
fully disordered initial state. At criticality, we have con- 
firmed that the correlation and response function behave 
as: C(t,t w ) = t- b f£{t/t w ) and R{t,t w ) » t" 1 " f%(t/t w ), 
with a = b = 2(3 /vz c . This conclusion has been sup- 
ported by a finite-size analysis on systems of linear sizes 
L = 10,20,30,40 and 50. From the datas on the corre- 
lation function we extracted the value X c /z c ~ 1.34 for 
the nonequilibrium critical exponent and checked its com- 
patibility with the response function datas. One may no- 
tice that this value is very close to that obtained in the 
3d-Ising model, that is A c ~ 2.8.0 It is remarkable to 
notice that the local scale-invariance prediction for the 
response function is very nicely fulfilled either at criti- 
cality and below. Although at criticality field-theoretical 
approach 14 showed a deviation from the local scale in- 
variance prediction,^] it is very difficult to see numer- 
ically this very small discrepancy. In the ordered phase, 
since the nonequilibrium exponent is very close to the 
free-field value d/2, we have shown that the autocorre- 
lation function has a scaling form which is very close to 
that of the free-field: C(t,t w ) = Mj q y~ x / Z with the scal- 
ing variable y = (x + \) 2 /{Ax) where x — t/t w . The 
direct calculation of the FDR X(t,t w ) showed that, as 
expected, in the low temperature phase it vanishes like 
Xoo(t,t w ) — t~ a fx(t/t w ) where fx is a scaling function 
simply related to the local scale-invariance scaling func- 
tions fc and /r. The same analysis at the critical point 
gave a limit FDR X^ = 0.43, which is in agreement, 
within the error bars, with the field theoretical value ob- 
tained in ref.^3] To conclude, we have studied the 3d— XY 
model in the context of aging and showed that the results 
obtained are consistent with the general picture usually 
drawn. Moreover, we have given a new verification of the 
local scale-invariance predictions, the first for a non-scalar 
order parameter. 
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